Superconductivity and crystalline electric field eff'ects in the filled skutterudite series 

Pr(Osi_,,Ruj4Sbi2 
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X-ray powder diffraction, magnetic susceptibility x{T), and electrical resistivity p{T) measure- 
ments were made on single crystals of the filled skutterudite series Pr(Osi_3;Rui,)4Sbi2. One end of 
the series (a; = 0) is a heavy fermion superconductor with a superconducting critical temperature 
Tc — 1.85 K, while the other end {x = 1) is a conventional superconductor with Tc ~ 1 K. The 
lattice constant a decreases approximately linearly with increasing Ru concentration x. As Ru (Os) 
is substituted for Os (Ru), Tc decreases nearly linearly with substituent concentration and exhibits a 
minimum with a value of Tc = 0.75 K at a; = 0.6, suggesting that the two types of superconductivity 
compete with one another. Crystalline electric field (CEF) effects in Xdc(T) and p(T) due to the 
splitting of the Pr^+ nine-fold degenerate Hund's rule J = 4 multiplet are observed throughout the 
series, with the splitting between the ground state and the first excited state increasing monotoni- 
cally as x increases. The fits to the Xdc(r) and p{T) data are consistent with a Fa doublet ground 
state for all values of x, although reasonable fits can be obtained for a Fi ground state for x values 
near the end member compounds (a; = or a; = 1). 

PACS numbers: 71.27.-|-a, 74.25.Fy, 74.25. Ha, 74.62. Dh 



I. INTRODUCTION 



The filled skutterudite compound PrOs4Sbi2 was re- 
cently discovered to be the first Pr-based heavy fermion 
superconductor, with a superconducting transition tem- 
perature Tc = F85 K and an effective mass m* « 50 mg, 
where nie is the free electron massji*^ Features in the 
dc magnetic susceptibility Xdc{T), specific heat C{T), 
electrical resistivity p{T), and inelastic neutron scatter- 
ing (INS) can be associated with the thermally depen- 
dent population of the ninefold degenerate Pr'^+ J = 4 
Hund's rule multiplet split by a cubic crystalline electric 
field (CEF). These data suggest that the ground state of 
PrOs4Sbi2 is a F3 doublet, separated from a Fs triplet 
first excited state by ~ 10 KpLSi The possibility of a Fi 
singlet ground state has also been put forward based 
on other measurements some of which also consider 
tetrahedral symmetry operators in their calculations of 
the CEF Hamiltonian of PrOs4Sbi2.^ It has been pro- 
posed that the superconductivity in PrOs4Sbi2 may be 
due to quadrupolar fluctuations,^ a claim that has been 
supported by /iSR^ and Sb-NQRi measurements, which 
indicate a strong-coupling isotropic energy gap of 2A « 
bksTc. Other intriguing effects are seen in PrOs4Sbi2, 
including multiple superconducting transitions^iSiiS and 
phases, and an ordered phase that is observed in high 
magnetic fields and low temperaturesii^ This high field 
ordered phase (HFOP), which is seen in measurements of 
p(r)^ C(r)^ magnetization M{T),^^ thermal expan- 
sion a(T)^ and magnetostriction X{T)^^ in a magnetic 
field H, as well as measurements of p{H) isotherms, ^'^ 
appears to be related to the crossing of the CEF energy 
levels in magnetic fieldsi^ii^ In addition, neutron diffrac- 
tion experiments* indicate the presence of quadrupolar 
effects in the HFOP, analagous to those seen in PrPh^i^ 



PrOs4Sbi2 has proven to be a unique compound, and will 
continue to provide a fertile area of research for many 
years. 



The isostructural compound PrRu4Sbi2 displays su- 
perconductivity below Tc « 1.0 K and possesses an elec- 
tronic specific heat coefficient 7 approximately 5 to 10 
times smaller than PrOs4Sbi2, identifying it at as a con- 
ventional metal, or at most a borderline heavy fermion 
metallic It was previously reported, based on measure- 
ments of Xdc{T), to possess a Fi ground state and a 
F4 triplet first excited state « 70 K above the ground 
state. A later measurement of p(T) also supported this 
CEF level scheme»i2, PrRu4Sbi2 appears to be a BCS-like 
weak-coupling superconductor, with an isotropic s-wave 
energy gap of 2 A « S/csTc, as determined from Sb-NQR 
measurements^ At the present time, no quadrupolar ef- 
fects or features resembling the HFOP seen in PrOs4Sbi2 
have been reported in PrRu4Sbi2. 



The substitution of PrRu4Sbi2 into PrOs4Sbi2 to form 
Pr(Osi_2;Rua;)4Sbi2 was undertaken to investigate the 
evolution of the superconductivity, the CEF energy level 
scheme, and the heavy fermion state with Ru-doping, 
and to investigate the relationship, if any, between these 
three phenomena. The present study focuses on mea- 
surements of xi^i'^) ^-nd p{x,T), which have revealed 
the x-dependencies of Tc and the splitting between the 
CEF ground state and the first excited state. We are 
also in the process of investigating the heavy fermion 
state via measurements of C(T), and the upper critical 
field Hc2{T) through measurements of p{T,H) (which 
will also reveal the a;-dependence of the HFOP), and will 
report these results in a future publication. 
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II. EXPERIMENTAL DETAILS 

Single crystals of Pr(Osi_a;Ru2;)4Sbi2 were grown us- 
ing an Sb flux method. The elements (Ames 5N Pr, Colo- 
nial Metals 3.5N Os and 3N Ru, and Alfa Aesar 6N Sb) 
were sealed under 150 Torr Ar in a carbon-coated quartz 
tube in the ratio 1 : 4-4a; : 4x : 20, heated to 1050 °C at 
50 °C/hr, then cooled at 2 °C/hr to 700 °C. The samples 
were then removed from the furnace and the excess Sb 
was spun off in a centrifuge. The crystals were removed 
from the leftover flux by etching with dilute Aqua Regia 
(HCliHNOaiHaO = 1:1:3). 

X-ray powder diffraction measurements were made at 
room temperature using a Rigaku D/MAX B x-ray ma- 
chine. The only significant impurities in any of the sam- 
ples were identified with free Sb that was still attached 
to the crystals. Each Pr(Osi_a;Rua;)4Sbi2 sample crystal- 
lized in the LaFe4Pi2 structure^" with a lattice constant 
a that decreased roughly linearly with increasing Ru con- 
centration X, as displayed in Fig. ^ A silicon standard 
was used in order to achieve a more accurate determina- 
tion of the lattice constant. Measurements of Xdc vs tem- 
perature T were made in a magnetic field H of 0.5 tesla 
between 1.8 and 300 K in a commercial Quantum Design 
superconducting quantum interference device (SQUID) 
magnetometer. Measurements of p and Xac were made 
as a function of T down to 1.2 K in a ^He cryostat and, 
for several of the samples, down to 0.1 K in a '^He-'*He 
dilution refrigerator. 

III. RESULTS 
A. Magnetic Susceptibility 

Displayed in the main portion of Fig. [21 is a plot of the 
dc magnetic susceptibility Xdc as a function of temper- 
ature T for single crystals of Pr(0si_2,Ru2,)4Sbi2 with 
various values of x. Above T « 100 K, the inverse mag- 
netic susceptibility 1/xdc is linear, indicating Curie- Weiss 
behavior. The data have been corrected for excess Sb by 
assuming that the high temperature effective moment, 
fics, of Pr should be equal to the Hund's rule free ion 
value of 3.58 /ib for Pr^+, where /is is the Bohr magne- 
ton. Any deviation from this value was attributed to 
free Sb inclusions in the Pr(Osi_a;Rua;)4Sbi2 crystals. 
The most significant effect on Xdc from this correction 
was not the small diamagnetic Sb signal but instead the 
change in overall scaling due to the difference in mass 
used to calculate Xdc in units of cm'^/mol from the raw 
magnetization data. The calculated percentages of mass 
attributed to Sb out of the total sample volume for all 
values of x are listed in Table The estimated value of 
the Sb mass depends slightly on the CEF ground state 
used to make the fit correction; only the values for a Fa 
ground state are given for simplicity. 

All of the Pr(Osi_a:Rua;)4Sbi2 samples exhibit features 
(peaks or plateaus) in Xdc that can be attributed to CEF 



effects. These features are the focus of the two insets in 
Fig. [3 The low temperature Xdc data for the samples 
from a; = to a; = 0.4 are shown in Fig. |2Ia), while Fig. 
|2Ib) similarly displays data for the samples from x — 0.5 
to X = 1. An explanation of the fits used to determine 
the CEF parameters from the Xdc data, as well as the 
parameters themselves, is given in section FlVBI 

Low temperature (< 2 K) ac magnetic susceptibility 
Xac vs T data for Pr(Osi_2;Rua;)4Sbi2 are shown in Fig. 13 
A sharp diamagnetic transition can be seen for all values 
of X, indicating the presence of superconductivity. The 
superconducting critical temperature Tc for each concen- 
tration was determined from the data displayed in Fig. |21 
as the midpoint of the diamagnetic transition. A plot of 
Tc vs X is displayed in Fig. [SJ and is discussed further in 
section llV Al An additional feature of note is the step- 
like structure that appears in the Xac data for PrOs4Sbi2 . 
Since double superconducting transitions have been ob- 
served in specific heat and thermal expansion measure- 
ments on both collections of single crystals and individual 
single crystalS)Si2ii& it is reasonable to assume that this 
step in the diamagnetic transition for PrOs4Sbi2 is also 
due to an intrinsic second superconducting phase instead 
of a variation of Tc throughout the multiple crystals used 
in the Xac measurements. None of the other concentra- 
tions display significant structure in their superconduct- 
ing transitions, although the transitions for x = 0.3 and 
X — 0.4 are much wider than for the other concentrations. 
This may be due to a variation of Tc between individual 
crystals for these two concentrations. 



B. Electrical Resistivity 

Fig. 21 displays high-temperature electrical resistivity p 
vs T data for Pr(Osi_xRu2;)4Sbi2 for various values of x 
between and 1. The values of p at room temperature, 
p{300 K), and the extrapolated values of p at zero tem- 
perature from fits to the p(T) data based on calculations 
of p{T) that incorporate CEF splitting of the Pr'^^ J = 4 
multiplet (see section llVB|) . p{0 K), are listed in Table 
13 Also listed in Table U is the residual resistivity ra- 
tio (RRR), defined as p(300 K)/p(0 K). It is surprising 
that the RRR of PrRu4Sbi2 is so much lower than that 
of PrOs4Sbi2, since they are both stoichiometric com- 
pounds and would be expected to have a low residual 
resistivity. A previous measurement of PrRu4Sbi2 found 
p(300 K) — 632 pil cm and an RRR of 25,^^ in reason- 
able agreement with the data presented in this paper. 
The low RRR of PrRu4Sbi2 is not presently understood. 

The electrical resistivity of Pr(Osi_:rRuj;)4Sbi2 below 
T = 2 K is shown in Fig. |S1 The data have been nor- 
malized to their values at 2 K in order to emphasize the 
superconducting transitions. The x — 0.7 sample did not 
display the onset of superconductivity down to the low- 
est measured temperatures and no data for this sample 
are shown in this plot; the heating due to large contact 
resistances in the x = 0.7 and x = 0.9 samples precluded 
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measurements below 1 K. The superconducting transi- 
tions as determined from p{T) are in reasonable agree- 
ment with those measured inductively (Fig. OJ, and the 
plot of Tc vs X is discussed in the following section. 



IV. ANALYSIS AND DISCUSSION 
A. Superconductivity 

The dependence of the superconducting transi- 
tion temperature Tc on Ru concentration x for 
Pr(0si_2:Rua;)4Sbi2 is shown in Fig. |B| Several concen- 
trations have more than one data point associated with 
them, which are from measurements of different crys- 
tals. These additional measurements were not shown in 
Figs. 0] and |S1 or listed in Table in the interest of clar- 
ity. The RRRs were nearly identical for all crystals of a 
given concentration, with the exception of the x = 0.2 
samples where the crystal with the lowest Tc in Fig. |H| 
had an RRR about half of that measured for the other 
two X = 0.2 samples, one of which is listed in Table 
The vertical bars in Fig. El £^re a measure of the width 
of the superconducting transitions, taken to be the 10% 
and 90% values of the resistance change associated with 
the transition. 

The trend of the Tc vs x data shown in Fig.^lsuggests 
a competition between the two different types of super- 
conductivity seen in PrOs4Sbi2 and PrRu4Sbi2. This 
competition suppresses Tc from both ends, culminating 
in a minimum of Tc = 0.75 K near x = 0.6. Specific 
heat measurements are in progress, and it will be inter- 
esting to see if the heavy fermion state can be correlated 
with Tc- The persistence of superconductivity through- 
out the series is unusual, as for heavy fermion f-electron 
superconductors both magnetic and nonmagnetic impu- 
rities generally produce relatively rapid depressions of Tc- 
When the impurity is of an element that would produce 
an isostructural superconducting compound, the trend is 
not as clear. For example, the Ui_2:La:rPd2Al3 system 
is similar to the Pr(0si_2:Ru2:)4Sbi2 system in that one 
end member compound, UPd2Al3, is a heavy fermion su- 
perconductor, while the other end member compound, 
LaPd2Al3, is a conventional BCS superconductor. Un- 
like Pr(Osi_a;Ru2;)4Sbi2, however, superconductivity is 
destroyed upon substitution on either end of the series. 
This persistence of superconductivity throughout of the 
Pr(0si_2:Rua;)4Sbi2 system for all values of x is ob- 
served in the CeCoi-ajIrj^Ins series of compounds, which 
is also superconducting for all values of a;.^^ This system's 
similarities to Pr(Osi_2;Rua;)4Sbi2 end there, because 
both end member compounds (CeCoIns and Celrlns) are 
heavy fermion superconductors in which the supercon- 
ductivity is believed to be magnetically mediated and to 
possess nodes in the energy gap A(k)ii2^ 

This nodal energy gap structure may be in contrast 
with PrOs4Sbi2, where ^SR^ and Sb-NQR'' measure- 
ments indicate an isotropic energy gap, a condition which 



could occur if the superconductivity in PrOs4Sbi2 was 
mediated by quadrupolar fluctuations. It is also gen- 
erally the case that superconductors with isotropic or 
nearly isotropic energy gaps are relatively insensitive to 
the presence of nonmagnetic impurities. Thus, the grad- 
ual decrease of Tc, and the presence of superconductivity 
for all values of x in Pr(Osi_xRua;)4Sbi2, provides further 
evidence for an isotropic energy gap and quadrupolar su- 
perconductivity in PrOs4Sbi2, since PrRu4Sbi2 also pos- 
sesses an isotropic superconducting energy gapiiS The 
minimum in Tc near a; = 0.6 could then be attributed to 
a shift from quadrupolar mediated heavy fermion super- 
conductivity to phonon mediated BCS superconductiv- 
ity. On the other hand, thermal conductivity measure- 
ments of PrOs4Sbi2 in a magnetic field have been inter- 
preted in terms of two distinct superconducting phases 
in the iJ — T plane, one with two point nodes in A(k) in 
low fields, and another with six point nodes in A(k) at 
higher fields. Since no thermal conductivity measure- 
ments were reported for fields below 0.3 T, the structure 
of A(k) is not known below this field. It is conceivable 
that, just as a magnetic field induces a change from a 
state with two point nodes into a state with six point 
nodes, the state with two point nodes is itself induced 
from an isotropic zero-field energy gap. Further mea- 
surements of the energy gap symmetry in zero and low 
magnetic field could shed light on this mystery. 



B. Crystalline Electric Field Effects 

The Xdc(T) and p{T) data for Pr(Osi_:rRux)4Sbi2 
were fit to equations including CEF effects, in a manner 
identical to that reported previouslyii^ The CEF equa- 
tions were derived from the Hamiltonian of Lea, Leask 
and Wolf (LLW).24 In the LLW formahsm, the CEF en- 
ergy levels are given in terms of the parameters xllw 
and W, where xllw is the ratio of the fourth and sixth 
order terms of the angular momentum operators and W 
is an overall energy scale factor. It was assumed that the 
CEF parameter y which controls the tetrahedral T^ crys- 
talline symmetry contribution to the Hamiltonian^ was 
small; thus, the calculations were made for a cubic Oh 
crystalline symmetry. Assuming that y is small implies 
that the main contribution to the crystalline electric field 
comes from the simple cubic transition metal sublattice 
(Os or Ru), as opposed to the more complicated tetrahe- 
dral Sb sublattice. The Xdc{T) data for x < 0.15 could 
be reasonably fit with either a Fs or a Fi ground state 
and a F5 first excited state. As x increases, the magni- 
tude of the peak in Xdc decreases more rapidly than the 
temperature T^ax at which the peak occurs. The peak 
also broadens until it resembles a hump. These changes 
with X make it unreasonable to fit a Fi-Fs CEF energy 
level scheme to the higher x data, since for these data 
an energy level scheme with the correct Tmax makes the 
peak too sharp, while the correct hump shape results in a 
Tmax that is too high. Thus, for the Pr(0si_2:Ru2.)4Sbi2 
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samples with x > 0.2, a ground state best approxi- 
mated the data. An example of a fit with a ground 
state for x — 0.6 is shown in Fig. Efa). A plot of the 
splitting between the ground state and the first excited 
state vs X is shown in Fig. |S1 including all reasonable fits 
of the Xdc{T) data. 

The Pr(Osi_j:Ru2;)4Sbi2 samples with x > 0.6 all dis- 
play upturns in Xdc{T) at the lowest temperatures (inset 
(b) in Fig.jSJ. If these upturns are due to the splitting of 
the CEF energy levels in a small magnetic field H, then 
it is expected that they would be more visible in the sam- 
ples with large x (more Ru than Os), where Xdc is smaller 
at low temperatures compared to the small x (more Os 
than Ru) data. The samples with x > 0.75, including 
PrRu4Sbi2, also display structure in these upturns that 
appear to be an additional peak near 5 K superimposed 
on the broad CEF hump, near the temperature of the 
CEF peak in PrOs4Sbi2 . The smooth progression of both 
the lattice parameter a and Tc indicates that there is 
no macroscopic phase separation of Pr(0si_2:Rua;)4Sbi2 
into PrOs4Sbi2 and PrRu4Sbi2. However, it is possible 
that the peak-like structure could be due to inhomoge- 
neous alloying of Os and Ru on an atomic scale, wherein 
each Pr'^^ ion sees a distribution of Os or Ru atoms, lead- 
ing to a variation in the CEF throughout the crystal. 
Unfortunately, this possibility would be difficult to es- 
tablish in the current experiments. The low-temperature 
upturn, especially in PrRu4Sbi2, could be attributed to 
either CEF splitting in H or paramagnetic impurities, 
both of which could produce a low-temperature increase 
in Xdc- 

Takeda et al. reported that PrRu4Sbi2 had a Pi sin- 
glet ground state and a r4 triplet first excited state, a 
CEF configuration that exhibits a plateau in Xdc at low 
temperaturesJii In the current experiment, the x = 0.9 
and X = 1 samples are the only ones in which a plateau 
is observed. In addition, while the other samples with 
X > 0.85 have their peaks reasonably well described by a 
Ps ground state, the fit predicts a saturation at T = K 
that is much lower than is observed in the data. However, 
the low-T upturn could be responsible for disguising both 
the maximum in x = 0.9 and x = 1 and the low tempera- 
ture saturation observed in the other high Ru concentra- 
tion samples. Accordingly, all the Pr(0si_2:Ru2;)4Sbi2 
data with x > 0.85 were fit assuming both a P3 — P5 
CEF energy level scheme and a Pi — P4 scheme, ignoring 
the low-temperature upturn; the x = 0.85 fits are shown 
in Fig.[7fb). Both fits are represented in the splitting be- 
tween the ground state and first excited state A_Egs_ics 
vs X plot of Fig. |S1 the results from all fits are also listed 
in Table n 

The electrical resistivity p{T) of Pr(Osi_i:Rux)4Sbi2 
was fit by a combination of scattering from impuri- 
ties, the atomic lattice (phonons), and temperature- 
dependent energy level populations due to the CEFpi^ 
The phonon contribution was represented by the mea- 
sured Plat of LaOs4Sbi2, an isostructural reference com- 
pound without f-electrons, for all values of x. This pro- 



cedure was validated by reproducing the results of Abe 
et al.^* with LaOs4Sbi2 instead of LaRu4Sbi2; as ex- 
pected, the Plat data of the two compounds appear to 
be nearly identical. The CEF contribution to p(T) con- 
sists of two terms, representing magnetic exchange and 
aspherical Coulomb scattering, which were assumed to 
be equally important when fitting the data;-^ Just as it 
was possible to fit p{T) of PrOs4Sbi2 with either a P3 
or a Pi ground state, all of the Pr(Osi_a;Rua:)4Sbi2 data 
were indifferent to the choice of either ground state. The 
splitting between the ground state and the first excited 
state (always a P5 triplet) was also nearly identical for 
fits with either ground state for a particular value of x. 
In the interest of simplicity, for the p{T) data, only the 
splitting between P3 and P5, is shown in Fig.|Sl 

The fit used to calculate Ai?3_5 for x = 0.15 is shown in 
Fig- 13b). It is evident that p{H) measurements will be 
required to elucidate the CEF ground state from trans- 
port measurements 

It is unclear what effect the CEF ground state may 
have on the superconductivity in Pr(Osi_a;Ru2:)4Sbi2. 
From a physical point of view, it is reasonable that an 
abrupt change in the ground state would produce an 
equally abrupt change in the physical properties. How- 
ever, it is difficult to conceive of a mechanism for this 
occurrence in the context of the LLW theory, since it is 
based on the interaction of the atomic lattice with a rare 
earth ion. If there is not an abrupt change in the lattice 
structure, one should not expect an abrupt change in the 
CEF ground state. It is therefore far more reasonable to 
consider a constant ground state, with the excited state 
varying as the Ru substitution changes the spacing of the 
atoms in the skutterudite lattice. The present data are 
most consistent with a constant P3 ground state, with 
the exception of the x = 0.9 and x = 1 data. However, 
when X is in the region 0.2 < x < 0.75 a P3 - P5 CEF 
energy level scheme is the only one which reasonably fits 
the XdciT,x) data. On the other hand, the possibility 
cannot be ruled out that this deep in a substituted sys- 
tem, a CEF analysis in the tradition of LLW may be 
unreliable due to the distribution of the two substituents 
(Or and Ru) in the near neighbor environment of each 
Pr'^^ ion. The gradual metamorphosis of the Xdc data 
does suggest that the CEF parameters are also changing 
gradually, but this may be misleading. Further experi- 
ments as well as theoretical analysis will be necessary to 
completely reveal the CEF ground state and its relation- 
ship to the superconductivity. 



V. SUMMARY 

The superconducting critical temperature Tc and crys- 
talline electric field (CEF) parameters of single crystals 
of Pr(Osi_a;Rua;)4Sbi2 have been deduced through mea- 
surements of x(T) and p(T) for < a; < 1 . The su- 
perconductivity, which is present for all values of x, ex- 
hibits a change in the sign of the slope in Tc{x) near 
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X = 0.6. The CEF ground state may also change from a 
ground state to a Ti ground state near this concen- 
tration, although more measurements are necessary to 
confirm this possibility. It is possible that the crossover 
from heavy fermion superconductivity that may be me- 
diated by quadrupolar interactions to nearly BCS super- 
conductivity occurs at this 'pseudocritical' concentration 

Xr,c = 0.6. 
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FIG. 1: Lattice parameter a as a function of Ru concentration 
X. The solid line is a linear least squares fit to a vs a;. 




FIG. 2: dc magnetic susceptibility Xdc as a function of tem- 
perature T between 1.8 and 300 K for single crystals of 
Fr(Osi_,rRua;)4Sbi2. Insct (a): Xdc vs T between 1.8 and 30 
K, showing the evolution of the peak due to crystalline elec- 
tric field effects for the Ru concentrations a: = to a: = 0.4. 
The X = 0.2 sample has been removed for clarity. Inset (b): 
as inset (a), but for Ru concentrations x = 0.5 to a; = 1. The 
samples with x = 0.75 and x = 0.95 were removed for clarity. 
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FIG. 3: ac magnetic susceptibility Xac as a function of 
temperature T between 0.1 and 2 K for single crystals of 
Pr(Osi_a;Rua;)4Sbi2. The data have been normalized to 
at T = 2 K and to -1 at T = K for clarity. 
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FIG. 4: Electrical resistivity p as a function of temper- 
ature T between 0.4 and 300 K for single crystals of 
Pr(Osi_a;Rua;)4Sbi2 with various values of x between and 
1. The samples with x = 0.05, x = 0.15, and x = 0.9 were 
removed for clarity. 
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FIG. 5: Electrical resistivity p as a function of temperature T 
between 0.4 and 2 K for single crystals of Pr(Osi-a;Rua,)4Sbi2 
with various values of x between and 1, normalized to their 
values at 2 K. The data for the sample with x = 0.7 is not 
shown because it did not superconduct down to the lowest 
temperature measured (see text for details). Similarly, the 
superconducting transition for x = 0.9 is not complete due to 
the limits of the experiment. 




FIG. 6: Superconducting critical temperature Tc vs Ru con- 
centration X for Pr(Osi-xRua;)4Sbi2. Filled circles: Tc ex- 
tracted from electrical resistivity p. Open squares: Tc deter- 
mined from ac magnetic susceptibility Xac- The straight lines 
are guides to the eye. 
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FIG. 7: Examples of CEF fits to the data, (a): dc mag- 
netic susceptibility Xdc(T) for x = 0.6 and x = 0.85 for 
Pr(Osi-a:Rua:)4Sbi2. The solid lines are fits assuming a Fa 
doublet ground state and a F5 triplet first excited state, and 
the dashed line is a fit assuming a Fi singlet ground state and 
a F4 triplet first excited state, (b): electrical resistivity p(T) 
for X = 0.15 between 1 K and 300 K. The fit is for a F3 ground 
state and a F5 first excited state. Fits with a Fi ground state 
were qualitatively identical, and so are not shown (see text for 
details). Inset to (b): p{T) for x = 0.15 between 1 K and 50 
K, displaying the low-temperature curvature in greater detail. 
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FIG. 8: The splitting between the ground state and 
first excited state Ai5gs-ies vs Ru concentration x for 
Pr(Osi_a;Rua:)4Sbi2, Calculated from fits of CEF equations 
to Xdc(r) and p{T), as described in the text. The solid line is 
a linear fit to Ai?gs-ies for a Fa doublet ground state and a 
Fs triplet first excited state calculated from the Xdc{T) data. 
For X < 0.15, a CEF energy level scheme with a Fi singlet 
ground state and a F5 first excited state also provided a rea- 
sonable fit to the Xdc(T) data, while a Fi ground state with 
a F4 triplet first excited state was also a possible energy level 
scheme for x > 0.85. 



TABLE I: Physical properties of Pr(Osi_i:Rua;)4Sbi2 com- 
pounds. X is the concentration of Ru; p(300 K) is the electri- 
cal resistivity p at 300 K; p(0 K) is p at K extrapolated from 
CEF fits (see text); RRR is the residual resistivity ratio, de- 
fined as p(300 K)/p(0 K); %Sb is the percentage of the mass 
attributed to free Sb in Xdc{T) assuming a F3 ground state; 
XIjUw and W are the Lea, Leask and Wolf parameters^— and 
AEa-b is the energy difference between ground state Fa and 
first excited state F;,. 
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